The ss and KK nature of / (980) in D s decays 

A. Deandrea( a \ R. Catto^, G.Nardulli^, A.D. Polosa^ and N. A. Tornqvist^ 
(a) Institut de Physique Nucleaire, Universite Lyon I, 43, bd du 11 Novembre 1918, F-69622 Villeurbanne Cedex, France 
(b) Departement de Physique Theorique, Universite de Geneve, 24 quai E.-Ansermet, CH-1211 Geneve 4, Switzerland 
(c) Dipartimento di Fisica, Universita di Bari and INFN Bari, via Amendola 173, 1-70126 Bari, Italy 
(d) Physics Department, POB 9, FIN -00014, University of Helsinki, Finland 

(December, 2000) 

We examine the D 3 — > /b(980)7T amplitude through a constituent quark-meson model, incor- 
porating heavy quark and chiral symmetries, finding a good agreement with the recent E791 data 
analysis of D s —* 3tt via /o(980). The /o(980) resonance is considered at the moment of production 
as an ss state, later evolving to a superposition of mainly ss and KK. The analysis is also extended 
to the more frequent process D 3 — * 4>ir. 

Pacs numbers: 13. 25. Ft, 12.39.Hg, 14.40. Cs 
LYCEN-2000-132 
UGVA-DPT/2000-11-1090 
HIP-2000-65/TH 



(i) Introduction 

The scalar mesons have remained a controversial enigma since a long time. There is today no consensus as to the 
true nature of especially the lightest scalars, the a, /o(980) and eto(980). Are these qq 0, tttt j|, KK bound state 
by hyperfine interaction || or multi-quark Q states? Is the expected scalar gluonium state |5||| present among these 
mesons, either as the dominant component, or through a small mixing? 

These are fundamental questions of great importance in particle physics. The mesons with vacuum quantum 
numbers are known to be crucial for a full understanding of the symmetry breaking mechanisms in QCD, and 
presumably also for confinement. The light cr, near 500 McV in mass and with a large width of the same order 
of magnitude, has reappeared on the scene (it has been on the short list of PDG since 1996 0J§]), and a growing 
number of analyses, using more sophisticated theoretical techniques, now find this elusive meson in 7T7t — > tttt, where 
its presence is almost hidden because of the Adler zero (in production experiments the a appears more clearly). 
Recently, in June 2000, a conference, entirely devoted to the cr, was held in Kyoto ||, and many new analyses of old 
data were presented, showing a wide agreement on a a pole near m — iT/2 — 500 — i250 MeV. 

Experiments studying charm decay to light hadrons are opening up a new unique experimental window for un- 
derstanding light meson spectroscopy, and especially the controversial scalar mesons. The scalars are abundantly 
produced in these decays, and in the recent E791 experiment jl0| the a contributes 46% of the D — > 3tt, and /o(980) 
over 50% of the D s -> 3tt Dalitz plot gj). 

D s decays are particularly interesting, since with the Cabibbo favored c — > s transition and the dominant spectator 
mechanism one expects the final state to be dominated by ss states. This mechanism is supported by the fact that, in 
D s decay, the well known nearly-ss vector state, the 0(1020), is abundantly produced, and one consequently expects 
that also the related scalar ss state should appear strongly. Since the /o(980) is produced copiously in D s decay, this 
supports the picture of a large ss component in its wave function. 

Previously this predominant ss nature of the / (980) has been supported by the radiative decay — * fo (980)7 
|L2| and by unitarized quark models although being just below the KK threshold it should also have a large 
component of virtual KK in its wave function. 

In this paper we shall assume that both the 0(1020) and the /o(980) are predominantly ss states, or at least that 
when they are produced in D s decay, the production is via the ss component. After production the ss core of /o(980) 
induces a virtual cloud of KK (due to the fact that /o(980) is just below the KK threshold and it couples strongly to 
KK) behaving as a large standing S— wave surrounding the relatively small qq core. We consider a model for heavy- 
light meson decays, the Constituent Quark Meson Model (CQM) [Q, so far successfully applied, and see if it predicts 
the decay rates for D s — > </>7r and /o(980)7r compatible with the recent data. The CQM has been widely exploited 
for heavy-light meson decays, in which the light quarks are u or d. Recently some of us Q studied D — > cr7r — > 3tt 
with the CQM, finding good agreement with the E791 data for this process and assuming the a is predominantly 
(uu + dd) I y/2. The role of the cr in B decays to 3 pions was also investigated |l5| . 

In a recent paper by Anisovich et ai, the authors adopt the hypothesis of an /o(980) having ss and (uu + dd)/^/2 
components with a mixing angle of —48° deduced from a phenomenological analysis of 0(1020) — * 7/o(980), 777, 777', 771- 
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and /o(980) — > 77 radiative decays (ll|. The iTA' component of /o(980) is neglected in this approach. On the other 
hand it has been shown by Markushin ^] how a KK molecular picture of /o(980) can explain the /o(980) — > rnr 
decay with no need of uu, dd components. As we shall see, our analysis favors a /o(980) produced as an ss state which 
evolves generating a virtual cloud of KK eventually decaying in an OZI allowed way to mr, as shown in Fig. 1. In 
fact we shall test the hypothesis that the ss component is substantial. 

(a) d s -> / 7r 

The CQM model can be extended in the strange quark sector solving the gap equation discussed in [ fL8| with a non 
zero current mass for the strange quark: 



II(m) = m - m - 8mGh(m 2 ) = 0, (1) 

id mo is tb 

proper time regularization: 



where G = 5.25 GeV"^ and toq is the current mass of the strange quark. The I\ integral is calculated using the 



- j^£_ I™ 9 rf4fc _ N c m 2 t to 2 to 2 \ 

The choice of the UV cutoff is dictated by the scale of chiral symmetry breaking A x = 4^/^ [[l9| and we adopt 
A = 1.25 GeV. The IR cutoff /i and the constituent mass to must be fixed taking into account that CQM does not 
incorporate confinement. This means that we have to enforce the kinematical condition to produce free constituent 
quarks M > toq + to, where M is the mass of the heavy meson and toq is the constituent mass of the heavy quark 
there contained. Considering that the heavy meson momentum is P M = wiqv^ + k^, being the heavy quark 
4- velocity and k^ the so called residual momentum due to the interactions of the heavy quark with the light degrees 
of freedom at the scale of Aqcd, the above condition coincides with v ■ k > to (since P = Mv, the 4-velocity of the 
meson is almost entirely carried by the heavy quark), or equivalently, in the rest frame of the meson, inf(fc) = to, 
meaning that the smallest residual momenta that can run in the CQM loop amplitudes are of the same size of the 
light constituent mass. The IR cutoff fj, is therefore (i ~ to. 

A reasonable constituent quark mass for the strange quark is certainly to = 510 MeV, considering the (f> meson as 
a pure ss state [0. Taking fj, = 0.51 GeV as an infrared cutoff, a value of to = 131 MeV, see Fig. 2, is required by 
the gap equation (consistently with the spread of values for the current s quark mass quoted into S). Varying the 
current strange mass in the range 60 — 170 MeV reflects into a small excursion of the constituent strange mass around 
the 500 MeV value. 

The free parameter of CQM is Ah defined by Ajj = Mh — tuq. The subscript H refers to the H— multiplet of 



Heavy-Quark-Effective-Theory (HQET) |20| H = (0~, 1~). In a similar way a As is associated to the S multiplet, 
S = (0 + , 1 + ). The latter is determined fixing Ah The related Ah, Ag values in the strange sector are shown 

in Table I. We consider the range of values Ah = 0.5,0.6,0.7 GeV in all numerical computations. This range is 
consistent with the condition M > toq + m. 

As a first test we compute the decay constant fo 3 given in CQM by the following expression: 

F 2^ZH{h + {AH + m)h{A H )) 

JD 3 = , — 1= , \i) 



where ttlh — to,£> s = 1.968 GeV and: 

iN c r° s d 4 l 



h(A) 



16 7T 3 / 2 7 1/A 2 S 3 / 2 

The renormalization constant Zh is given by: 



16tt 4 J {I 2 - m 2 )(v ■ I + A + ie) 

e -s(,n 2 -A 2 ) (l +er f(Av^)). (4) 



N r /•' "'' ds ....2 ,2 



Numerically we find: 



Z H * = (A H + m)^^- + h(A H ). (5) 



f Ds = 297115 MeV, (6) 
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(the error is computed varying Ah in the range of values quoted above) that is in good agreement with the value in 
the PDG §: 



Let us now consider the D & 

(h{q h )\A^ c) {q)\D s {p)) 



f D , = 280 ± 19 ± 28 ± 34 MeV. 
/o(980) semileptonic amplitude: 



(J) 



(p+q fo r 



Ju 



l fo 



''D, 



F (q 2 ), 



(8) 



with Fi(0) = i*b(0). This amplitude can be represented by the diagram in Fig. 3. CQM allows to model the /o 
vertex, indicated with a black spot, through the diagrams in Figs. 4 and 5 respectively. The former gives what we 
call the polar contribution to the form factors: considering the 0~ intermediate polar state one can compute Fq while 
1+ is connected to F\. The direct diagram, depicted in Fig. 5, gives access to both the computation of F\ and F . 
The method and the computation technique have been fully explained in |14|. Actually, in a factorization scheme, 
the amplitude describing the decay D s — > /o(980)7r is expressed by the product of two matrix elements. One is the 
semileptonic matrix element (H), the other is the well known: 



(*\Afa(q)\VAC)=iUg». 



(9) 



The g M in this matrix element selects Fq as the relevant form factor for the computation of the amplitude. The polar 



and direct contributions to the form factors must be added in order to determine Fo(q 



0). In the case at 



hand we use a coupling of /o(980) to the light quarks that, for SU3— flavor symmetry, is gf oSS = V2g aqq — 2.49\/2 
pl[ . This of course comes from the hypothesis of an /o (980) having a ss structure. If on the other hand one would 
adopt the picture given in [^6| , gf oss should be reduced by a factor of sin(</>), with 4> — —48°. The numerical value 
for Fq is obtained through the CQM expressions containing the 7i(A#) and ij(As) integrals discussed in JlJ]. The 
result is the following: 



F (g 2 = 0) = F (pol) 
The expression for the decay amplitude is: 



(0) 



G 



^ (dir) (0H0.64±r o i 



9D afo « = (f 7r+\H cS \Df) = -^V*V udai F (0)(m 2 Da -mjJU, 



(10) 



(11) 



where H c h is the effective Hamiltonian of Bauer, Stech and Wirbel [£2|, with a\ = 1.10 ± 0.05 fitted for D decays 
|E3|. Its numerical value is: 



g Dsfo7r = (2.08 



+ ait) x 10" 6 GeV, 



and therefore the predicted width is: 

T(Df - / (980)^+) = 



167rr7i|j 



^\(ml s ,m%,ml) = (3.27t° f 5 ) x 10" 



GeV, 



(A is the triangular function), to be compared with the PDG Q one: 

T(Df -> /o(980)vr+) = (2.39 ± 1.06) x 10" 14 GeV. 



(12) 



(13) 



(14) 



Considering the mixing of ss with the (uu + dd)/\/2 component, i.e., using the reduced coupling g/ oSS \/2sin(</i) 
discussed above, one would obtain a different prediction for the width: 



T{D S / (980)7t) = (1.86±8;H) x 10 



-14 



(15) 



Even if (^3|) and (|l5|) seem both to agree with the experimental value (|l4|), for reasons that will be explained soon, 
the CQM model definitely favors the ss scenario. 
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In order to make a comparison with E791 results |Tlj], we compute the branching ratio B(Df — > /o(980)-7r + — > 37r) 
estimating the coupling gf g7r+jr - for /o — > 7r + 7r~ through the following formula: 



47T /™2 



(16) 



where C = 2/3 x 0.68 ~ 4/9 is the fraction of the total /o(98Q) width related to the process /o — > 7r + 7r~ (the iT_fir 
decay mode is also considered, 2/3 being the isospin factor and 0.68 = r(/ — > 7T7r)/(r(/ — > 7T7r) +r(/ — > KK)) gj) 
and To = 44 MeV is the central value of the width of /o(980) found by E791 |pj| . Using the expression: 

r W , (980>„* 3.) = I * r^.) x ; ( ._^T4 ( ^w 

where the T's are computed assuming a pure ss component in production as in ( |l3| ) but substituting — ► s, and 
the co-moving width is: 



T f0 (s)=r ^4^L. (18) 



Assuming the E791 value m/ = 975 MeV |Tlj and taking the central value of the PDG branching ratio B(D S — > 
37r) = 1.0% (disregarding the large given uncertainty ±0.4) one finds that the fraction of D s decaying into 37r via 
/o(980) is given by: 

B(D+ -> /o(980)7r+ -> 3vr) = 50±j? %. (19) 

This turns out to be in good agreement with the E791 results Jil"| : 

B(D+ -► / (980)tt + -> 3tt) = (56.5 ± 4.3 ± 4.7) %, (20) 

as far as the central value of the experimental branching ratio B(D S — ► 37r) is concerned. Unfortunately the lack of 
precision in this measurement weakens the comparison between the CQM calculation and the E791 results. 

(Hi) D s — > (fnt 

In order to overcome the ambiguity between what found in the pure ss hypothesis and in the case of an /o(980) 
pictured as a mixture of ss and (uu + dd)/\^2 (both results (13) and (|l5| ) seem to be consistent with the experimental 
value (|l4|)) we consider the CQM analysis of the D s — > (fin decay in order to compare the widths T(D S — > <pTr) and 
T(D S — > /o7r). We believe that the ratio T(D S — > foir)/T(D s — > 07r) is the most reliable theoretical output of our 
approach to be compared with experimental data since it is less dependent on CQM parameters and experimental 
normalizations. We will find that also the T(D S — ► 4>tt) turns out to be larger with respect to the experimental 
one, analogously to what found in ( |l3| ) with respect to (|lj). Interestingly the ratio of the two widths is in very 
good agreement with the experimental ratio. This agreement is instead destroyed if a large mixing with uu and dd 
components is taken into account. 

Considering the 0(1020) as an ss state, we can again make use of factorization hypothesis p2| to compute the 
D s — > 4>tt amplitude. The semileptonic form factors are defined by: 

(cb(e,p')\A^ c) (q)\Df(p)) = (mns+m^Attfy 



t {£ *\ q) d P + p'y i A2{ q 2 )-(e*- q ) 2 -^^ ' — ' ~ 
(m Ds + TTifp) q A 



-(p + pTMq 2 )- (e*-q)^q^A 3 (q 2 )-A (q 2 )). (21) 



To avoid the singularity in q 2 = 0, the condition: 

A (0)=A 3 (0), (22) 

must hold. Moreover: 

Mq 2 ) = ^±^±A^) - ^^±Mq-). (23) 
2rri0 2rri0 
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Observing that q = (p — p') and using eq. (||) one obtains: 

9d.*k = (<f>Tr + \H cfi \Df) = ^V^VudaMml^m^e* ■ q)f„. (24) 

In order to compute the width, a sum over final <j> polarizations is performed. The form factor Aq(Q) can be computed 
with CQM following the approach outlined in pi| j. There is a direct and a polar contribution to the form factors, as 
in the diagrams in Figs. 4,5. The 4> particle attached to the loop is introduced, via Vector-Meson-Dominance (VMD), 
through an interpolating current: 

2 

J» = -^f 7m, (25) 



where the leptonic decay constant f$ is defined by: 

(VAC|V^(e)> = (26) 

with — — S7 M S, and: 

= -pcos^c), (27) 
V6 

= 39.4° being the 4> — uj mixing angle, see e.g. p5| . This current appears in the CQM loop integral calculation 
in correspondence of the vertex (light quark)- {(f)- (light quark). The lepton decay constant can be extracted from 
the value of the width: 

/ \ 2 

4na 2 



r(0 - e+e") = — \^fj OT0> (28) 

(as is checked with a straightforward calculation), given in the PDG We will use the value = 249 MeV 2 . 
Following p4| , the direct contribution to A o (0) is given by: 

fil ( m^LO — ) + 171^2 + 

V m D s J 



At\q 2 = 0) = -^yfZ^ e 

U c <t> 



rriD.m,!, 



where: 



7 I 2 r l . 

Z I m — m h mmsuj 

rriD 3 



(29) 



in 



ri = 2 ' 

and the integrals flj, which are functions of Ai, A2, x and w, are tabulated in |24f| and computed using x = m^, 
the w given above and the Ah values from Table I. A2 are here substituted by Ah — m^uj. The polar form factors 
needed to compute ^4q 01 (<7 2 = 0) are, according to eqs. ( p2| ) and ( |23| ) , A^ ol (0) and A^iO), given respectively by [pi| : 

A? ol (Q) = ^f g ^ + (C - 2^) (30) 

= -v2, 9 ^F+^ ( r Ds+m0) 

where now a) = Tp^S <JV = -f^- and F + is defined by: 



(31) 

m D s (l+) 
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(VAC|A M |D a (l+)) = i^m Ds{1+) v^F + , (32) 
and computed in analogy with (^), see Jig] . The £ and \i strong coupling constants are described in H): 

■ y/z^(-Cl 1 -2— +mz) (33) 



V2gvf4>C4> V m 

V2m2 



C = t 1 - V^jjZs (m fi 2 + 2fi 3 + ^4 + 5 - m 2 Z) , (34) 

where the Qj integrals involved are functions of Ai = Ah, A 2 = As, x — m<p and u) — (Ai — A 2 ) /m^. The 
renormalization constant Zg is defined by: 

Z- 1 = (A s - m)^|^ + J 3 (A S ). (35) 

Of course everywhere m = 510 MeV is the constituent mass of the strange quark, /i and £ are affected by a considerable 
uncertainty due to varying Ah in the range of Table 1; this is reflected in a large uncertainty in the polar form factor 
ylg Ol (0) = — 0.16lg 2 7 . The direct form factor is instead more stable against Ah variations being ^4q"'(0) = 1.26^o'}5. 
Considering that: 

A (0) = ^ ol (0) + ^n' (0 ), (36) 
we can readily compute the CQM ratio of widths in the case of only ss in the production: 

to be compared with the PDG one || : 

R = r( ^t — = 0-49 ± 0-20- (38) 

r(u+ -+ 4m+) v ' 

(iv) Conclusions 

If one adopts the hypothesis of Anisovich et ai, the ratio R would be reduced to R = 0.22 ± 0.12. E791 also 
measures the ratio T(D S -> 37r)/r(D s -> / (980)tt) = 0.245 [0 with a very small uncertainty. If one considers ( |l7j ) in 
the limit of narrow width for /o(980), one obtains T(D S — > /o(980)7r — * 37r) = C r(D s — > /o(980)7r) (see the discussion 



after eqjl6|). On the other hand E791 finds that T(D S -> / (980)tt -> 3tt) = 56.5%r(D s -> 3?r) or, in other words, 
they measure R = 0.62 with a very small error. This indeed agrees with the known PDG result (pi]). 

The computed widths £> s — > /o7r and Z? s — > 4>n are both larger with respect to the corresponding experimental 
values, nevertheless their ratio is only 20% smaller than the experimentally estimated ratio. Our results favor the 
scenario of an /o(980) made of an ss core surrounded by a standing S'-wave of virtual KK. A large uu, dd component 
in /o(980) seems also excluded by the fact that in D — > 3n decays the /o(980) is weakly produced poj . 

Therefore, in conclusion, our work supports a description of Jo (980) as ss, with virtual KK cloud. Any substantial 
mixture of uu, dd seems excluded. Light-quark phenomenology has been fighting since a long time to understand the 
scalar mesons. It now appears that heavy meson decays may be able to clarify this difficult problem: nemo propheta 
in patria. 
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K 

Fig. 1 - After the production at short distances in D s 
decay via its ss component (see Figs. 3-5), the /o(980), 
being just below the KK threshold, evolves in time gen- 
erating a substantial KK component (with larger spatial 
dimension than ss) that can decay, in OZI allowed way, 
to 2tt. 
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Fig. 2 - Gap equation zero for a fixed value of mo- The 
masses are expressed in GeV. Here mo = 131 MeV is the 
current mass of the strange quark. 




fo 

Fig. 3 - The semileptonic amplitude. The vertex with 
the weak current and the fo can be modeled with CQM 
as is described in Figs. 3 and 4. The same diagrams with 
4> in place of fo are also considered. The 4> resonance is 
introduced, via VMD, through an interpolating current 




fo 



Fig. 4 - The polar diagram. The polar contribution to 
the form factor is reliable when computed near the pole 
mass. The uncertainty in the extrapolation q 2 — > re- 
flects in the violation of the condition Fq° (0) = Ff° (0). 
This kind of uncertainty is taken into account in our cal- 
culation. 



l+k 




fo 

Fig. 5 - The direct diagram. The condition ^ o dir (0) = 
^, dil '(0), avoiding the spurious singularity in q 2 = in 
(g), is automatically satisfied. 
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